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We report on theoretical studies of recently discovered microwave-induced resistance oscillations 
and zero resistance states in Hall bars with two occupied subbands. In the same results, resistance 
presents a peculiar shape which appears to have a built-in interference effect not observed before. 
We apply the microwave-driven electron orbit model, which implies a radiation-driven oscillation of 
the two-dimensional electron system. Thus, we calculate different intra and inter-subband electron 
scattering rates and times that are revealing as different microwave-driven oscillations frequencies 
for the two electronic subbands. Through scattering, these subband-dependent oscillation motions 
interfere giving rise to a striking resistance profile. We also study the dependence of irradiated 
magnetoresistance with power and temperature. Calculated results are in good agreement with 
experiments. 

PACS numbers: 



I. INTRODUCTION 

Transport excited by radiation in a two-dimensional 
electron system (2DES) is currently a central topic from 
experimental and theoretical standpoints*. The interest 
is focussed not only on the basic explanation of a physical 
effect but also on its potential applications. In the last 
decade it was discovered that when a Hall bar (a 2DES 
with a uniform and perpendicular magnetic field (B)) is 
irradiated with microwaves, some unexpected effects are 
revealed, deserving special attention from the condensed 
matter community: microwave-induced (MW) resistance 
oscillations (MIRO) and zero resistance states (ZRS)2r— . 
These remarkable effects show up at low B and high mo- 
bility samples, specially ZRS where ultraclean samples 
are needed. Different theories have been proposed to ex- 
plain these striking effects^— but the physical origin is 
still being questioned. To shed some light on the physics 
behind them, a great effort has been made, specially from 
the experimental side, growing better samples, adding 
new features and different probes to the basic experimen- 
tal setup, Of course the experimental results 
always mean a real challenge for the existent theoreti- 
cal models. Thus, a comparison of experiment with the- 
ory could help to identify the importance of the invoked- 
mechanisms in these theories. 

One of the most interesting setups, carried out recently, 
consists in using samples with two or three occupied 
subbands 23 . These samples are either based in a dou- 
ble quantum well structure or just one single but wide 
quantum well. The main difference in the longitudinal 
magnetoresistance (R X x) of a two-subband sample is the 
presence of magneto-intersubband oscillations (MISO) 28 . 
These oscillations occur due to periodic modulation of 
the probability of transitions through elastic scattering 
between Landau levels (LL) of different subband. The 
MISO peaks corresponds to the subband alignment con- 
dition A = nhw Cl where A is the subband separation and 



w c the cyclotron frequency. Because of elastic scattering 
of electrons between LL of different subband, the proba- 
bility rate is maximal under this condition. Under MW 
irradiation the first experimental results 2 ^ of R xx showed 
the interference of MISO and MIRO without reaching the 
ZRS regime. Later on, further experiments realized at 
higher MW intensities and mobility samples, showed the 
MW-response to evolve into zero resistance states for the 
first time in a two occupied subband sampled. In the 
same experiment 2 ^ it was also observed a peculiar R xx 
profile with different features, regarding the one-subband 
case 2 - - — , affecting only valleys and peaks of MIRO's in a 
surprising regular way. Thus, in valleys we observe two 
nearly symmetric shoulders, one at each side of mini- 
mum which could correspond to a more intense transport 
through the sample. On the other hand, in the peaks we 
observe narrower profiles, regarding again the one sub- 
band case, meaning a smaller transport. 

In this article, we theoretically study magnetoresis- 
tance of a Hall bar being illuminated with MW radiation 
when two electronic subbands participate in the trans- 
port. We apply the theory developed by the authors, the 
MW-driven electron orbits mode^^—, which we extend 
to a two-subband scenario. According to this theory, 
when a Hall bar is illuminated, the electron orbit cen- 
ters of the Landau states perform a classical trajectory 
consisting in a harmonic motion along the direction of 
the current. Thus, the 2DES moves periodically at the 
MW frequency altering dramatically the scattering con- 
ditions and giving rise eventually to MIRO and ZRS. In 
some cases the transport is reinforced producing MW- 
induced R xx peaks; in others, transport is weakened giv- 
ing rise to valleys. In a double subband scenario the sit- 
uation gets more complicated but with a richer physics. 
On the one hand, due to the presence of MW, we have 
two 2DES (twoo subbands) moving harmonically at the 
MW-frequency. On the other hand, we have two possible 
scattering processes with charged impurities: intra and 
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II. THEORETICAL MODEL 




FIG. 1: Schematic diagram for the wide quantum well and 
the corresponding electronic wave functions. ^1 and are 
the individual Fang-Howard wave functions for the left and 
right triangular-like built-in potential wells, "fs and ^Pyt are 
the symmetric and antisymmetric wave functions of the wide 
quantum well. The quantum well has a width of 45nm as in 
the experiments. 



The MW driven electron orbits model, was developed 
to explain the R xx response of an irradiated 2DEG at 
low B. We first obtain an exact expression of the elec- 
tronic wave vector for a 2DES in a perpendicular B, a 
DC electric field and MW radiation which is considered 
semi-classically. Then, the total hamiltonian H can be 
written as: 



H = 



^ + \m*w 2 c {x - X) 2 - eE dc X 
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X is the center of the orbit for the electron spiral motion: 
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Eq the intensity for the MW field and Ed c is the DC 
electric field in the x direction. H\ is the hamiltonian 
corresponding to a forced harmonic oscillator whose orbit 
is centered at X. Hi can be solved exactl y 33 ' 34 , and using 
this result allows an exact solution for the electronic wave 
function of H to be obtained^ 3 ^— : 



^> N (x,t) OC (f> n (x - X - Xd(t),t) 



(3) 



inter-subband. We then calculate the two correspond- 
ing elastic impurity scattering rates, obtaining that the 
intra is, approximately, three times larger than the in- 
ter. This means first, that the current is mainly sup- 
ported by intra-subband scattering processes. Secondly 
and more important, the competition between intra and 
inter-subband scattering events under the presence of ra- 
diation alters significantly the transport properties of the 
sample. This is reflected in the R xx profile through a 
strong and peculiar interference effect. As in experi- 
ments, our calculated results recover the presence of new 
features regularly spaced through the whole MIRO's pro- 
file, mainly two shoulders at minima and narrower peaks. 
We identify such features with situations where the in- 
terference is constructive and the current is reinforced 
(shoulders around minima) meanwhile in other cases the 
interference is destructive giving rise to a less intense 
current (thinner peaks). Within the same theory, we 
have obtained also ZRS in the same position of exper- 
iments and with the same MW-frequency dependence. 
Finally, we have studied the influence of MW-frequency 
(w), MW-power (P) and temperature (T) on MIRO's of 
the two subband sample and the obtained results are also 
in reasonable agreement with experiment^,. 



where <p n is the solution for the Schrodinger equation 
of the unforced quantum harmonic oscillator, x c i(t) is 
the classical solution of a forced and damped harmonic 
oscillator 
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where 7 is a phenomenologically-introduced damping fac- 
tor for the electronic interaction with acoustic phonons. 

Then, the obtained wave function is the same as the 
standard harmonic oscillator where the center is dis- 
placed by x c i(t). Thus, the electron orbit centers are 
not fixed, but they oscillate harmonically at w. This 
radiation — driven behavior will affect dramatically the 
charged impurity scattering and eventually the conduc- 
tivity. Thus, we introduce the scattering suffered by the 
electrons due to charged impurities. If the scattering is 
weak, we can apply time dependent first order perturba- 
tion theory. First, we calculate the impurity scattering 
rat o 5 ' 32 ' 36 between two oscillating Landau states \&jVj fi- 
nal state 1 i' m (x, t): 
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FIG. 2: Schematic diagrams of electronic transport corre- 
sponding to peak and valley scenarios and intra and inter- 
subband scattering types. When the MW field is on, the or- 
bits are not fixed but oscillate at w. In 2a, in the intra-band 
processes corresponding to peaks, the electronic orbits are go- 
ing backward giving a larger advanced distance and current. 
Yet, only 3/4 of these processes develop through this channel. 
In 2.b we represent the inter-subband scattering processes for 
a peak, where due to the lower scattering time the orbits are 
now going forwards, yielding a smaller advanced distance and 
current. According to our calculations they are 1/3 of the 
total. Summing up all process, the total advance distance is 
given by |AX^T a + jAX," 1 ^. Then, we obtain less current 
with respect to the one-subband case due to the lower contri- 
butions of the inter-subband events. This result correspond 
to the obtained narrower profiles at maxima. The valleys 
situation reflected in 2c and 2d, can be explained in similar 
terms as the peaks but now the discussion and results are 
going to be the opposite. Then, considering jointly all scat- 
tering processes, intra and inter-subband, we obtain, at both 
sides of minima, more current than the one-subband case due 
to larger contribution of the inter scattering events. This is 
the physical origin of the two shoulder that can be observed 
a both sides at minima. 



trie constant, and qo is the Thomas- Fermi screening 
constant^ 7 .. 

After some lengthy algebra we arrive at the expres- 
sions for the intra-subband W* n ?I a and the inter-subband 



scattering rates: 
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where rii is the density of impurities, T the width of 
the Landau states, A12 the subband separation and Fs 
and Fa are the form factors. To obtain the form factor 
expressions we have considered, as in experiments^ 3 -, a 
highly doped wide quantum well with ~ 9.1 x 10 n cra. 
In this type of wells, as more electrons are added, their 
electrostatic repulsion forces them to pile up near the well 
sides and the resulting electron charge distribution ap- 
pears increasingly as bilayer (see Fig. 1). In other words, 
and effective electrostatic barrier is built up in the middle 
of the well separating the two electron layers in GaAS. In 
the case of a double quantum dot the barrier is made of 
ALGaAs or AlAs. As a result we obtain at each side of 
the wide quantum well a potential profile similar to the 
inversion layer of a 2DES that, in a good approximation, 
can be considered triangular close to the heterojunctions 
(see Fig. 1). Next, we have to obtain first the corre- 
sponding wave functions of these triangular-like built-in 
potential profiles. Then, due to its great simplicity and as 
a first approach, we have applied the Fang- Howard varia- 
tional treatment (see reb 37 ' 38 ) that proposes as electronic 
wave function (Fang-Howard wave function): 
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where V s is the scattering potential for charged 
impurities 3 - 7 -: 

y. = E 9<? /, (6) 
^ 2Se{q + q a ) 

S being the surface of the sample, e the GaAs dielec- 



where b is the corresponding Fang- Howard variational pa- 
rameter. According to this simple but efficient approach, 
b results to be mainly dependent on the two-dimensional 
charged impurity densit y 37 ' 38 . Now starting from this 
variational wave function we can build ^s(A) which are 
the corresponding symmetric (antisymmetric) wave func- 
tion of the wide quantum well(see Fig. 1). Finally the 
form factors are obtained: 
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FIG. 3: a) Calculated R xx vs B for dark and MW situa- 
tions; the ZRS is marked with an arrow, b) Same as a) for 
2-subbands and 1-subband. We observe clearly the new fea- 
tures showing up in the 2-subband curve comparing to the 
1-subband; shoulders at minima and narrower peaks. In the 
inset the ZRS region is blown up. Shoulders and narrower 
peaks are the outcomes of the interference between the intra 
and inter-subband scattering processes. 



-q(z 
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where q is the electron wave vector exchanged in the scat- 
tering. We have supposed a symmetrical delta doping, 
being d the average separation between the impurities 
and the 2DES at each side of the wide quantum well, 
i.e., spacer distance. This distance depends on the sam- 



ple. In practical terms it varies between 10nm to even 
larger than lOOnm. In our calculations we have used a 
numerical value of d = 70— 90nm. According to the Fang- 
Howad variational approac h ' 3S , where the parameter b 
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is given by: 



33m*ni 
8h 2 e 



1/3 



(12) 



, and applying the experimental sample parameters^,, 
we have calculated a numerical value for b ~ 0.25nm . 
Consequently the average thickness of the triangular 
wells (z) ~ 11 — 12nm, where (z) is related with b by 
(z) = |. Following again with the experimental param- 
eters at hand^ in terms of impurity density, well thick- 
ness, etc., we have made an averaged estimation of the 
relative values of F$ and Fa resulting in 



3.2 x \F A 



(13) 



where we have used an average value for q = 0.5nm _1 
that is of the order of the Fermi wave vector in agreement 
with the experimental impurity density. Substituting the 
obtained form factors in the scattering rates we can even- 
tually reach averaged values for those rates that result to 
be related by: 



(14) 



Here we have considered that the cosine average value, 
• for A 12 > hw c and we have carried out 



sg7rAij 
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->• 



the sum Y^T-\ - 

Once we know the intra and inter-subband scattering 
rates, we consider that when an electron undergoes a 
scattering process jumping from the initial state to the 
final one, it takes an average time 



i^~intra(inter) 



1 



t t 7 intra (inter) 
W n ,m 



(15) 



Following the model described in ref.— , we next find the 
average effective distance advanced by the electron in ev- 
ery scattering jump in the presence of radiation AX MW 
, generalizing the previous results to a two subbands sce- 
nario: 
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FIG. 4: Calculated R xx vs B for different frequencies. In 
a) for lower frequencies, from 50 to 100 GHZ, and in b) 
for higher, from 110 to 180 GHZ. We observe the pres- 
ence of ZRS in all curves with shifting positions depend- 
ing on the frequency and the interference features in peaks 
and valleys. ZRS positions and shifts are in agreement with 
experiment a 23 ' 40 and are similar to the one-subband case. 



av-MIC 

intr a(inter) 



AX° + Acos(w(T intra ( inter) )) (16) 



where AX° is the effective distance advanced when there 
is no MW field present. Applying the important previous 
result of 
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we can write the final expression for the total average 
distance advance due to both kinds of scattering, intra 



and inter, AX%%: 
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A COS [w(T intr a)] + A COS [w(T inter )] (19) 

A cos —(T in ter) + A cos [w(r mter )] (20) 



This significantly alters the scattering conditions regard- 
ing the one-subband case mainly affecting MIRO's peaks 
and valleys. 



Finally the contributions intra and inter-subband to 
the longitudinal conductivity can be calculated: a xx oc 

/ dE[(^^) intr a + (^^)i„ter](fi-ff), being fi and 
ft the corresponding distribution functions for the initial 
and final Landau states respectively and E energy. The 
obtained final expression for the conductivity is given by: 
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where Xs = 2 \^" bT and E\ and i?2 are the energies 
of the first and the second subband respectively. Finally, 
to obtain R xx we use the relation R xx — , , ~ 

where cr^ ~ ^ and cr^ < ct X!/ . 

The expression of the conductivity a xx shows the phys- 
ical equivalence to a situation with only one scattering 
time and two different oscillations frequencies for the 
MW-driven subbands: w/3 for the intra-band scattering 
and w for the inter. They demonstrate also the origin for 
the regular and strong interference profile observed in ex- 
periments where the factor 1/3 is essential to obtain the 
interference effect. A different factor would produce a to- 
tally distinct interference and also distinct R xx response. 
This factor comes from the calculation of the squared 
magnitude of the corresponding form factors Fi ntra and 
Fintra which eventually determine the different scattering 
rate between the intra-subband and the inter-one pro- 
cesses. These form factors depend mainly on the vari- 
ational parameter b and on the averaged wave vector q 
which subsequently are determined by two-dimensional 
impurity density rij. Therefore we can conclude that the 
crucial parameter 1/3 and eventually the obtained inter- 
ference profile will be mainly dependent on rij. During 
the scattering jump the electron perceives an approxi- 
mately three times faster MW-driven oscillation of the 
2DES when is inter-subband with respect to the intra- 
subband. This equation reflects also the important result 
that the intra-band conductivity (upper term in the a xx 
expression) is three times larger than the inter-subband 
(lower term of a xx ). Then, the total current is mainly 
supported by intra-subband scattering processes regard- 
ing the inter processes in a relation of approximately 
three to one. 

Based in these results, we can explain physically how 
the interference between both types of scattering process 
work, producing an excess of current at minima and a 
lack of current at maxima. In Fig. 2 we present schematic 
diagrams for the different situations. In all of them the 
MW field is on and the electronic orbits are not fixed, 
and instead move back and forth through x c i. In Fig. 
2a, the intra-subband scattering corresponds to a peak 
and the orbits moves backwards during the jump, then 
on average, electrons advance further than in the no MW 
case. Thus, we obtain more current giving rise to peaks 
because the average advanced distance is directly pro- 
portional to the conductivity. Yet, according to our cal- 
culations, only three out of four of the total scattering 



processes take place through the intra-subband channel. 
The other one out of four are inter-subband processes (see 
Fig. 2b). In them, due to the slower scattering time (or 
higher w) the scattering jump occurs when the electronic 
orbits are moving forward. Therefore, the average ad- 
vanced distance is smaller than in the intra processes. 
Eventually at the peaks, the total advanced distance, 
§AX*f t ^ + \AX^ r , is smaller than the one-subband 
case. Accordingly, this is directly translated to the ob- 
tained current, being reflected in the narrower peaks pro- 
file (see Fig. 3b). 

The valleys situation can be explained in similar terms 
as the peaks. Now, in the intra processes (see Fig. 2c) the 
scattering jump takes place when the electronic orbits are 
moving forward giving a smaller electronic advanced dis- 
tance and current. This situation gives rise to valleys and 
eventually if the MW power is big enough the ZRS can 
be achieved. As before, these intra-subband scattering 
processes correspond to three out of four of the total (see 
Fig. 2c). The other one out of four develops through the 
inter-subband channel (see Fig. 2d). In these processes 
, again due to the smaller scattering time or higher w, 
the scattering jump occurs when the electronic orbits are 
going backwards giving rise to a larger average advanced 
distance and current (see Fig. 2d). Then, summing up 
all scattering processes we obtain more current than the 
one-subband case. This is the physical origin of the two 
shoulder that can be observed a both sides at minima of 
R xx versus b. 



III. RESULTS 

In Fig. 3a, we present calculated R xx vs B for dark and 
MW situations and frequency / = w/2n = 143 GHz. We 
can observe MISO for the no-MW curve, MIRO for the 
MW curve and the ZRS marked with and arrow. In Fig. 
3b, we present the same as in 3a, but for 2-subbands and 
1-subband cases. The later has been obtained making 
A12 — > 0. Contrasting both curves we observe the new 
features appearing regularly spaced in peaks and valleys: 
two nearly symmetric shoulders in valleys and narrower 
peaks regarding the 1-subband curve. According to our 
model, these new features are result of the interference 
between the competing intra and inter-subband scatter- 
ing processes. This interference effect is mainly based 
in the different scattering rates between intra and inter- 
subband scattering events. On the one hand this is going 
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FIG. 5: a) Calculated results of power dependence of R xx vs B. MW power decreases from 4 mW to darkness and MIRO's 
decrease too. This is a similar behavior as the 1-subband result, b) Calculated A.R XX = R xx w — R xx vs power, for data coming 
from peaks (1) and (2) of the upper panel. The obtained fits for both peaks mean a sublinear P-dependence in agreement with 
previous experimental 22 and theoretical results. 



to be obviously reflected in the different scattering times: 
the inter process is three times smaller than the intra- 
subband. On the other hand, in the different capability 
to support the current. The intra-subband processes are 
able to support three times more current than the inter. 
Thus, in valleys, we observe a constructive interference 
effect giving rise to two shoulders meaning more current 
through the sample, meanwhile the narrower peaks mean 
a destructive interference and less current. The presence 
of ZRS is explained similarly as in ref£. 

In Fig. 4 we present calculated R xx vs B for different 
MW frequencies. In 4a, for a lower frequencies range 
(from 50 to 100 GHZ) 40 and in 4b, for a higher fre- 
quencies range (from 110 to 180 GHZ) 23 . We observe 
in all cases the presence of ZRS with shifting position 
depending on / and with reasonable agreement with ex- 
periment. All curves present the peaks and valley fea- 
tures meaning the importance of the interference effect 
which shows up independently of /. This /-dependence 
of ZRS positions was previously and similarly observed 



in one-subband samples 2 ^. 

In Fig. 5a we present P-dependence of R xx vs B for 
/ = 143 GHz. We observe that MIRO's decrease as P, 
(and Eq), gets smaller from 4 mW to darkness in similar- 
behavior as the 1-subband results^. In 5b, we present 
AR XX = Rx 4 x W ~ Rxx vs f° r data coming from peaks 
(1) and (2) of Fig. 5a, where R xx is the magnetoresis- 
tance for darkness and R xx w is the magnetoresistance 
when the MW field is on. We fit the data obtaining for 
both peaks a sublinear P-dependence, R xx oc P a where 
a < 1 and explained in terms of: 

E oc VP => Rxx oc VP (22) 

and in agreement with current 23 and previous 
experimental 2 ^ and theoretical 4 ^ results. 

In Fig 6a. we present the T dependence of R xx vs B 
for / = 143 GHz. As in 1-subband samples, we observe 
a clear decrease of MIRO for increasing T, eventually 
reaching a R xx response similar to darkness. In Fig. 6b, 
we present lnAR xx vs X/T for data coming from peaks 
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FIG. 6: a) Calculated results of T dependence of R xx vs B for / = 143 GHz. As in 1-subband samples, we observe a clear 
decrease of MIRO for increasing T, eventually reaching a R xx response similar to darkness, b) lnAR xx vs 1/T for data coming 
from peaks (1) and (2) of 6a. The fitted curves show the relation R xx cx T~ 2 in agreement with experiment. The T-dependence 
is explained with the damping parameter 7 which represents the interaction of electrons with acoustic phonons. 7 is linear 
with 232»22, thus an increasing T means an increasing 7 and smaller MIRO's. When the damping is strong enough (higher T), 
MIRO's collapse 



(1) and (2) of Fig. 6a. The T-dependence, according to 
the model, is explained with the damping parameter 7 
which represents the interaction of electrons with acous- 
tic phonons. 7 is linear with T— thus an increasing T 
means an increasing 7 and smaller MIRO's. When the 
damping is strong enough (higher T) MIRO's collapse. 
The curves of Fig. 6b show the relation R xx cx T~ 2 in 
agreement with experiment^. 



IV. CONCLUSIONS 

In summary, we have theoretically studied the recently 
discovered microwave- induced resistance oscillations and 
zero resistance states in Hall bars with two occupied 
subbands. MW-driven magnetoresistance presents 
a peculiar shape which appears to have a built-in 



interference effect not observed before. Applying the 
microwave-driven electron orbit model, we calculate 
different intra and inter-subband electron scattering 
rates under MW, revealing that the first is three times 
greater than the latter. This is physically equivalent 
to different microwave-driven oscillations frequencies 
for the two electronic subbands. Through scattering, 
these subband-dependent oscillation motions interfere 
giving rise to such a striking resistance profile. In the 
valleys the interference is constructive giving rise to two 
symmetric extra shoulders at each side of minima. In the 
peaks the interference is destructive producing smaller 
peaks regarding the one-subband case. The factor ^ 3 is 
essential to obtain this strong and regularly spaced in- 
terference effect. We study also the dependence on MW 
frequency, MW intensity and temperature. Calculated 
results are in good agreement with experiments. 
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